Two-gap model for underdoped cuprate superconductors 
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Various properties of underdoped superconducting cuprates, including the momentum-dependent 
pseudogap opening, indicate a behavior which is neither BCS nor Bose-Einstein condensation (BEC) 
like. To explain this issue we introduce a two-gap model. This model assumes an anisotropic 
pairing interaction among two kinds of fermions with small and large Fermi velocities representing 
the quasiparticles near the M and the nodal points of the Fermi surface respectively. We find that 
a gap forms near the M points resulting into incoherent pairing due to strong fluctuations. Instead 
the pairing near the nodal points sets in with phase coherence at lower temperature. By tuning the 
momentum-dependent interaction, the model allows for a continuous evolution from a pure BCS 
pairing (in the overdoped and optimally doped regime) to a mixed boson-fermion picture (in the 
strongly underdoped regime). 
PACS numbers:74.20.De, 74.20.Mn, 71.10.-w 



The underdoped cuprates are characterized by a pseu- 
dogap opening below a strong doping (S) dependent 
crossover temperature T* (SY, above the superconducting 
critical temperature T C (S) Q. By decreasing the doping 
the temperature T* increases, while the superconduct- 
ing critical temperature T c decreases until the insulating 
state is reached. The different behavior of T* and T c as 
doping is varied, finds a counterpart in the different be- 
havior of the coherence energy scale, obtained in Andreev 
reflection measurements and the single-particle gap, 
observed both in angle-resolved photoemission (ARPES) 
and in tunneling experiments This has triggered a 
very active debate on the relevance of a non-BCS super- 
conductivity and of a BCS-BEC crossover in these ma- 
terials |^-^|. In particular ARPES shows that below T* 
the gap opens around the M points of the Brillouin zone 
[i.e. (±7r, 0), (0, ±7r)] suggesting that T* can be inter- 
preted as a mean-field-like temperature where electrons 
start to form local pairs without phase coherence. How- 
ever, it is also found 0] that below T* substantial por- 
tions of the Fermi surface remain gapless. This behavior 
can be described neither by BCS nor by Bose-Einstein 
condensation (BEC) schemes. Instead it is suggestive of 
strong pairing between the states around the M points 
and of weak coupling near the zone diagonals. Various 
other experiments ||]9| carried below T c show a doping 
and temperature dependence of the gap anisotropy and 
therefore are again suggestive of a strong anisotropy in 
the pairing potential. 

In this letter we explore a new direction (neither BCS 
nor BEC) focusing on the consequences of a strongly 
anisotropic interaction. To this aim we introduce a two- 
gap model, where strongly paired fermionic states can 
coexist and interplay with weakly coupled pairs in dif- 
ferent regions of the Fermi surface (FS). This line of 
thinking was partly explored in Refs. |^,^[, where only 
the extreme strong-coupling limit of one component was 



considered. In particular the view of Ref. [pi would al- 
low to describe only the very underdoped region of the 
cuprate phase diagram. Our approach, instead, turns 
out to be sufficiently flexible to investigate with continu- 
ity the evolution of the bifurcation (taking place around 
optimal doping in the real systems) between a mean-field- 
like pairing temperature T* and the coherence supercon- 
ducting temperature T c . We implement our model by 
a two-band system with different intraband and inter- 
band pairing interactions. One band has a large Fermi 
velocity vf and a small attraction giving rise to largely 
overlapped Cooper pairs with weak superconducting fluc- 
tuations. On the contrary, the other band has a small vf 
and a large attraction resulting into tightly bound pairs 
having strong fluctuations. At variance from the mod- 
els of mixed fermions and bosons, we keep the fermionic 
nature of both the weakly and strongly bound Cooper 
pairs. 

A possible realization of a strongly momentum- 
dependent effective interaction in underdoped cuprates 
has been proposed in connection to the occurrence of a 
charge instability for stripe formation [ p~0|JTl] ] . It was in- 
deed suggested that the tendency to spatial charge order 
(which evolves into a spin-charge stripe phase by lower- 
ing the doping) gives rise to an instability line. This line 
T^ tripe (S) starts from a Quantum Critical Point (QCP) at 
T = near optimal doping and increses by lowering the 
doping. By approaching the instability line the pairing 
is mediated by the strong attractive quasi-critical stripe 
fluctuations, which affect the states on the FS in a quite 
anisotropic way: 
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where U is the residual repulsive interaction between the 
quasiparticles, 7 is a damping parameter, and q c is the 
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wavevector of the stripe instability. The crucial parame- 
ter n 2 is a mass term proportional to the inverse square 
of the correlation length of charge order £~ 2 and provides 
a measure of the distance from criticality. At T = 0, in 
the overdoped regime, k 2 is linear in the doping devia- 
tion from the critical concentration k 2 oc (8 — 8 C ). On 
the other hand, in the finite-temperature region above 
8 C n 2 oc T. In the underdoped regime k 2 vanishes ap- 
proaching the instability line T^ tripe {8) and extend the 
singular potential to finite temperatures. For n 2 k 0, 
the fermionic states around the M points are such that 
kp — kp' ~ q c and interact strongly. These are the so- 
called "hot spots" they have a low dispersion, and 
possibly form tightly bound local pairs giving rise to the 
pseudogaps below T* > T^ tripe (S). On the other hand, 
"cold" states in the arcs of FS around the zone diago- 
nals T-Y or T-X (nodal points) have larger dispersions 
and interact more weakly since V e g is now cut-off by q c . 
In the underdoped regime k 2 ~ at higher and higher 
temperatures by lowering the doping and V c g has a more 
dramatic effect. On the contrary, in going to the opti- 
mum and the overdoped region V e g is cut-off first by the 
temperature and then by the doping itself. All the states 
then interact more isotropically. 

- The two-gap model. Irrespectively of the origin of 
the anisotropy, in the presence of a strong momentum de- 
pendence both of (i) the effective pairing interaction and 
(m) the Fermi velocity, we must allow for enough freedom 
of the pairing and of its fluctuations in order to capture 
the relevant physical effects of the anisotropy. Follow- 
ing the above discussion we introduce a simple two-band 
model for the cuprates. We describe the quasiparticle 
arcs of FS about the nodal points by a free electron band 
(labelled below by the index 1) with a large Fermi veloc- 
ity vpi = kpi/rrii and the hot states about the M points 
with a second free electron band, displaced in momentum 
and in energy from the first, with a small vf 2 — fcra/^ 
(see Fig. 1) 0. 
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Fig. 1: Sketch of the Fermi surface of underdoped cuprates 
with quasiparticle arcs (thick solid line) and patches of quasi- 
localized states (thick dotted line) and the Fermi surface of 
the two-band spectrum (thin solid line). 

The assumed electronic spectrum consists therefore of 
two different free electron dispersions, £i(p) = p 2 /2m\ 
and £2(p) =| p Po | 2 /2rn2 + £o, displaced by a mo- 
mentum po and by an energy Eq ~ Ef\ introduced 



to allow the chemical potential to cross both bands: 
Epi — Eq + Ep2- To connect our two-band structure 
with the single-band dispersion of the cuprates, we choose 

po = (±7r,0),(0,±7r). 

This choice gives rise to two branches for the band 2 
along the x and y directions. Moreover, we only con- 
sider Cooper pairs of zero total momentum formed by 
time-reversed momentum eigenstates. Therefore the 2- 
2 pairs are always formed by (k, — k) states on portions 
of the band 2 symmetrically located on opposite sides 
with respect to the T = (0, 0) point. If the pairs have 
a s-wave symmetry, the branches along x and y of the 
band 2 are equivalent and just one can be considered. 
On the other hand, in the case of d-wave pairing, the 
pairs in the branch along x have a different phase from 
the pairs in the branch along y and both branches have 
to be treated. Since in this paper our main interest is the 
interplay between strongly and weakly coupled pairs ir- 
respectively of their symmetry, for simplicity we consider 
the s-wave problem . The model Hamiltonian for pairing 
in the two-band system is taken to be 

H =^2E kl n kal +^V lJ (k,k / )c+ +p1j ct k , lj C-kiiCk+pii 

hai kk'pij 

(2) 

where i and j run over the band indices 1 and 2 and a 
is the spin index. The interaction is approximated by a 
BCS-like attraction given by 

Viji^k') = - gij Q{uj - | fc(fc) |)e(a,o- I Ci(fc') I) (3) 

with an energy cutoff ujq. The strongly g-dependent effec- 
tive interaction in the particle-particle channel V^ff(q — 
k — k') of the original single-band system is accounted 
for by the 2x2 scattering matrix g. The matrix ele- 
ments gij couple the electrons within the same band (gn 
and 1722) and between different bands (gi 2 = g 2 i)- The 
self-consistency equation for the superconducting fluc- 
tuation propagator p[ in the matrix form is given by 
L = g + gtlL, where the particle-particle bubble opera- 
tor for the two-band spectrum has a diagonal 2x2 matrix 
form with elements Ilii(g) and 1122(g) V&j - The resulting 
fluctuation propagator is given by 

f (n) _ (hi -n u (g) gi2 V 1 . ( ] 

where we have defined g^ = {g^ r )ij- It turns out use- 
ful to define the temperatures Tjl and T° 2 as gn — 
n n (0, T) = g u - Pl ln^= Pl In g 22 - tt 22 (0, T) = 

g 22 - p 2 \n^ = p 2 hi^r, where pi = mi/(2ir) is the 
density of states of the i-th band. In the underdoped 
regime, to emulate the hot and cold points related for 
instance to the system near a stripe instability, we as- 
sume the following relations between the gij elements: 
522 ^ V/k, 2 » gn ~ V/\q c \ 2 ~ g 12 . Then one has 
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hi — 1/SU) 322 ^ I/522, 512 - -512/(311522)- In this 
limit and T c ° 2 (with T c ° 2 » T^) assume the value 
of the two BCS critical temperatures for the two decou- 
pled bands. The mean-field BCS superconducting critical 
temperature for the coupled system Tj? is defined by the 
equation detL _1 (q = 0,T C °) = 0. We obtain T c ° > T c ° 2 
given by 



T c ° = JTC>T c ° 2 exp 



45L 
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- The Ginzburg- Landau approach. The role of fluctu- 
ations can be investigated within a standard Ginzburg- 
Landau (GL) scheme, when both p2522^o < Ep2 and 
wo < Ep2- Under these conditions the chemical poten- 
tial is not affected significantly by pairing. Moreover, in 
order to remain within the GL approach, we will assume 
that fluctuations from the BCS result are not too strong. 
The relevance of the space fluctuations of the order pa- 
rameter is assessed by the gradient term coefficient 77, 
which provides the momentum dependence of the fluc- 
tuation propagator in Eq.(Q). In particular the expan- 
sion of the particle-particle bubbles, in terms of q, reads 
U n (q) - nn(0) - pimq 2 and U 22 (q) ^ n 22 (0) -P2V2Q 2 - 
Here r\i(i — 1,2) is the gradient term coefficient of the 
7-th band which, in 2D and for a free electron band, is 
given by m = {7({3)/32ir 2 )v 2 Fl /T 2 , with m > m @. 
In the absence of the interband coupling gi 2 , 771 provides 
the (large) gradient term coefficient corresponding to (the 
weak) superconducting fluctuations for the band with a 
large vfi, while (the small) r\2 corresponds to (strong) 
superconducting fluctuations for the band 2. For the cou- 
pled system near T®, the coefficient 77 in terms of 771 and 
r\2 is obtained by evaluating (L^ 1 )^ oc (e + r/q 2 ) in terms 
of the relative temperature deviation e = (T — T®)/T®. 
We get the expression 



V = 



Pi (522 - n 22 )?7i + p 2 (5n - 1111)772 

-T( 92 2 - n 22 )dnn/dr - T{g X i - nii)dn 22 /dT' 



(6) 



where all the bubbles are evaluated at q = 0. Using 
the definitions for 2 and the condition deti _1 (q = 
0, Xj?) = 0, the coefficient 77 can be explicitly written as 
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am + "2772, with — = 2 , t0;t0 v 

«2 pip 2 ln (TU/T^) 



(7) 



and ot\ + Q!2 = 1. The presence of a fraction of electrons 
with a large 771 increases the stiffness of the whole elec- 
tronic system (i.e. increases 77 with respect to 772). How- 
ever when the mean-field critical temperature T® is much 
larger than the correction to 772 due to the interband 
coupling is small. At the same time the Ginzburg num- 
ber is large implying a sizable mass correction Se(T) to 
the "mass" e(T) of the bare propagator L(q). The renor- 
malized critical temperature TJ, given by the equation 



e(T c r ) + (5e(T c r ) = 0, is lower than T c ° 0. By evaluating 
the renormalized gradient term coefficient rf in the pres- 
ence of the mass correction, we find that this is still given 
by Eq. (|) with T c ° replaced by T r c . Therefore rf = r}(T r c ) 
is greater than r](T®). This result indicates that the mass 
renormalizations of the fluctuation propagator tends to 
lower T c and, at the same time, increases the gradient 
term coefficient 77 by increasing the coupling to rji. As 
a consequence the effective Ginzburg number is reduced 
and the system is stabilized with respect to fluctuations 
allowing for a coherent superconducting phase even in 
the limit 772 — ^ 0. Within the GL approach we associate 
the temperature Tj? ~ T® 2 to the crossover temperature 
T* and XjT to the superconducting critical temperature 
T c of the whole system. In the region T c < T < T* 
the pseudogap is formed in band 2. Superconducting 
fluctuations only affect band 2 while are immaterial for 
band 1 , where the Fermi surface is mantained until phase 
coherence sets in. Within the Stripe-QCP scenario the 
coupling (722 is related to the singular part of the effec- 
tive interaction mediated by the stripe fluctuations. 522 
is the most doping dependent coupling and attains its 
largest value in the underdoped regime, where k 2 van- 
ishes at higher temperatures. The regular parts of the 
interaction 511 and 1712, being cut-off by q c , are instead 
only weakly doping-dependent. In the region of validity 
of the GL approach, the explicit calculations show that 
r{5) = (T c ° - T c r )/r c ° ~ (T* - T c )/T* is increasing by in- 
creasing (722, «.e., by decreasing doping. For small values 
of r we find that both T* and T c increase. This regime 
corresponds to the overdoped and optimally doped re- 
gion. Specifically, above optimum doping, g 2 2 and 511 
become comparable and the two lines merge together. 
For r ~ 0.25 -j- 0.5, T c is instead decreasing while T* is 
always increasing by decreasing doping. The large val- 
ues of r, which are attained in the underdoped region 
show that we are reaching the limit of validity of our 
GL approach. We think, however, that the behavior of 
the bifurcation between T* and T c represents correctly 
the physics of the pseudogap phase, while a quantitative 
description would require a more sophisticated approach 
like a RG analysis. 

- The strong- coupling limit. In the very low doping 
regime, where T* has increased strongly, the value of 
(722 can be so large to drive the system in a strong cou- 
pling regime for the fermions in band 2 (^2322^0 > Ep 2 ). 
In this case, taking ojq > Ep2 the chemical potential is 
pulled below the bottom of the band 2. In this limit 
of tightly bound 2-2 pairs, the propagator of the super- 
conducting fluctuations in band 2, i.e. 7^22(9), assumes 
the form of a single pole for a bosonic particle. Since 
Epi is still the largest energy scale in the problem, the 
fermionic character of the particles in band 1 is preserved. 
The critical temperature of the system is again obtained 
by the vanishing of det L^ 1 (q = 0) where, however, the 
chemical potential is now self-consistently evaluated in- 
cluding the selfenergy corrections to the Green function 
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in band 2 and the fermions left in band 1. 
n 22 = /O 2 wo/|^ 2 | and 



T c = T c ie xp 



IM2||Mb| 



PlP2^0 (|M2| - \Vb\ 



One gets 



(8) 



where /i 2 is the chemical potential measured with respect 
to the bottom of the band 2 and = p 2 <? 22 wo repre- 
sents the bound-state energy. The calculation of II 22 at 
small q leads to a finite value of r/ 2 = l/(8m 2 |/x 2 |), while 
the small-o expansion of det-L -1 provides the new r\ co- 
efficient [n8f 



V = Vi+ f^T In 2 (T^/T^)r l2 



(9) 



In this strong-coupling case most of the non-mean-field 
effect has been taken into account by the formation of the 
bound state occurring at a very high temperature of the 
order /? 2 <? 22 wo, which provides the new T* in this regime. 
r\ ~ 771 stays sizable and the fluctuations will not strongly 
further reduce T r c with respect to T c °: T c ~ T r c ~ T c °. In 
this low-doping regime T 7 ^ Tc approaches its largest val- 
ues before T c vanishes. 

The strong coupling limit of our model shares some 
similarities as well as some important differences with 
phenomenological models of interacting bosons and 
fermions [^J. In particular, in the model of Ref. || 
pairs of non interacting electrons are scattered from the 
Brillouin zone near the nodal points into dispersionless 
boson states, localized about the M points at an energy 
Eb- The correspondence with our model can be seen 
by noticing that the tightly bound bosonic states cor- 
respond to p2g22 UJ o 3> Ep 2 and gn = 0. The tightly 
bound dispersionless bosonic states are fully incoherent 
(i]2 = 0), while the fcrmionic states are unpaired as long 
as bosons and fermions are independent. The fermion- 
boson coupling gyi effectively introduces an e-e coupling 
of the order of pip 2 woffi 2 /£s and drives the system su- 
perconducting. In this particular limiting case of our 
model we recover the results of Ref. with an explicit 
expression for the bosonic level Eb = \p\ — while in 
Ref. Q it is a free phenomenological parameter. 

- Conclusions. In our paper we have analyzed the 
pairing properties of the underdoped cuprates in terms 
of an effective two-gap model, motivated by the strong 
anisotropy of the band dispersions and of the effec- 
tive pairing interaction. The crucial schematization was 
based on the introduction of two bands weakly coupled in 
order to preserve a substantial distinction between the su- 
perconducting order parameter in different regions of the 
momentum space. This has to be contrasted with more 
standard approaches producing one single gap (even with 
complicated momentum dependence) and a common fluc- 
tuating order parameter all over the FS. Our approach 
allows for different fluctuation regimes for pairs in dif- 
ferent k regions. According to our analysis, the strongly 
bound pairs forming at an high temperature T* can ex- 
perience large fluctuations until the system is stabilized 



by the coupling with less bound BCS-like states, leading 
to a coherent superconducting state at T c < T* . T c and 
T* merge around or above optimum doping, where g 22 , 
according to our stripe scenario, becomes of the order 
of <7n. Our model shares similarities with the fermion- 
boson models for cuprates [f|D to which it reduces in the 
strong-coupling limit for g 22 3> g±i ~ 0. The two-gap 
model considered here applies to a much wider doping 
region and is more suitable to describe the crossover to 
the optimal and over-doped regime, where no preformed 
bound states are present and the superconducting tran- 
sition is quite similar to a standard BCS transition. 
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